Transition temperature of ferromagnetic semiconductors: a dynamical mean field 

study 
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We formulate a theory of doped magnetic semiconductors such as Gai-ajMn^^As which have at- 
tracted recent attention for their possible use in spintronic applications. We solve the theory in the 
dynamical mean held approximation to find the magnetic transition temperature Tc as a function 
of magnetic coupling strength J and carrier density n. We find that Tc is determined by a subtle 
interplay between carrier density and magnetic coupling. 
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Diluted magnetic semiconductors have attracted much 
recent attention [Q for their potential use in spintronic 
devices. The prospect of carrying out both information 
processing and storage on the same chip is an excit- 
ing possibility. For applications it is desirable to find 
materials which are ferromagnetic at as high a temper- 
ature as possible, so the discovery ||] of a ferromag- 
netic transition with Tc as high as llOK in MBE grown 
Gai_a;Mn2:As, has inspired a great deal of interest, fn 
addition to its potential technological significance, the 
ferromagnetism of Gai_a;Mna;As is an important funda- 
mental condensed matter problem. The cause of ferro- 
magnetism in Gai_a;Mna;As (and similar materials, e.g 
Ini_^Mn2:As where Tc « 25 — SOifjEl) is controversial: 
the different proposed mechanisms ppH] do not qualita- 
tively agree with each other. In this paper we present 
a new theoretical approach which allows calculation of 
magnetic transition temperatures (and other properties) 
over a much wider temperature range than had previ- 
ously been possible and provides new insights into the 
factors controlling Tc- 

It is well established that in III-V systems such as 
Gai_xMn2,As, the Mn ions go in substitutionally at the 
cation (Ga) sites and contribute itinerant holes to the 
GaAs valence band. The experimental hole density n is 
typically a small fraction (10% or so) of the Mn concen- 
tration perhaps due to strong localization at As antisite 
defects so the Gai_xMn2:As system could be considered 
partially compensated. The Mn ion has a half filled d- 
shell and acts as a S* = 5/2 local moment; the itinerant 
carriers are locally magnetically coupled to the Mn spins 
via an exchange coupling J which may be either anti- 
ferromagnetic (J > 0) (for electrons) or ferromagnetic 
(J < 0) (for holes : physical Gai-^^Mn^jAs case). Q. 

It is generally accepted p]-p^ that the magnetic semi- 
conductors are described by a generalized Kondo lattice 
model: 
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where u labels the relevant bands of the semiconductor 
(the two hole bands in hole-doped GaAs for example), 
Vr is a potential arising from randomness in the host lat- 
tice (e.g. As antisite defects), Ri are the positions of the 
Mn dopants, W is the (presumably coulombic) potential 
arising from the Mn dopant, Jaf is a direct antiferro- 
magnetic exchange between Mn spins arising from other 
orbitals unrelated to the doped holes, J is the local ex- 
change coupling between the spin of the Mn and the spins 
of the semiconductor carriers. We normalize the S func- 
tion in the J term to the volume — 5.65A per GaAs 
unit. The large (S — 5/2) value of the Mn spin justi- 
fies treating the spins classically, so that the partition 
function Z may be determined by finding the free energy 
T({Si}) of holes in a fixed spin configuration and then 
averaging over spin configurations, i.e. 
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The key issue is therefore the evaluation of T'dSi}) . 
From Eq. 0, we see that F is the free energy of nonin- 
teracting carriers in a spin dependent potential, which 
may have randomness both from the distribution of Mn 
positions and from spin disorder. Further, as will be 
shown explicitly below, the relevant temperatures are 
small compared to the hole Fermi energy {Ef\ so that 
T is to a good approximation simply the carrier ground 
state energy in the given spin configuration. The crucial 
quantity governing T^ is the change in F as the spin con- 
figuration goes from disordered to ordered. As we will 
show below, the change in F involves several competing 
effects not evident in previous calculations of Tc such as 
the static mean field theory [|. 
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To obtain more detailed information, we consider here 
the idealized Kondo lattice model, in which Jaf , Vr and 
W in Eqj^ are neglected. These may easily be included 
in our formalism, and will be discussed in a future paper. 
We absorb the magnitude of the spin into the definition 
of J and study first the T ^ limit of the fully polarized 
ferromagnetic state. Si || z at T = 0. In this state the 
carriers feel a spatially varying spin dependent potential 
with mean strength xJ per GaAs unit cell. It leads to 
a shift in the band offset (upwards for one species and 
downwards for the other) proportional to xJ. For val- 
ues of J less than a critical value Jc, the spin dependent 
potential does not lead to any bound states. The wave 
functions are scattering states and the band offset is es- 
sentially ±xJ. The energy is given simply by filling the 
up and down bands up to the appropriate chemical po- 
tential. The critical value Jc corresponds to a magnetic 
coupling strong enough to bind a hole of the appropriate 
spin to a Mn site, causing a 'parallel-spin' impurity band 
to split off from the main band. Sanvito et al. Q used 
the local spin density approximation (LSDA) and super- 
cell methods to study Gai-a^Mn^^As with a dilute but 
spatially ordered Mn lattice. The LSDA prediction for 
J was found to be very close to the critical value needed 
for impurity band formation (Fig. 11 of Ref. jl^], shows a 
band shift, within 15% of AE — lOey/Mn, implying in 
our conventions J = leV); experimental estimates tend 
to be somewhat lower. 

We now consider energetics of arbitrary spin configura- 
tions. In the small J limit, straighforward perturbation 
theory shows that 

Ji^.= i^jMi?.,^.)S.-S, (3) 

with x{R±jRj) the static spin susceptibility computed 
from Eq |l| with J = 0. This is often referred to ||-||] 
as the 'RKKY' limit although strictly speaking the term 
RKKY refers to the behavior of x at distances long com- 
pared to the spacing between carriers and Eq. ^ applies 
even for spins closer together than this distance. When 
Eq. H applies, the ordering wavevector is the one which 
maximizes x and Tc ^ J^, which is also the static mean- 
field result Q. Note that for n greater than a (numeri- 
cally small) critical value ric, the maximum in % is at a 
nonzero wavevector, leading to a non-ferromagnetic or- 
dered state. 

In the J —^ oo limit, at all times each carrier is bound 
to an Mn site with a binding energy proportional to J 
and spin parallel to the Mn spin on that site. The de- 
pendence of energy on spin configurations arises because 
in the paramagnetic state some hopping processes are 
blocked |ll] and is therefore set by the impurity band- 
width which is never large (because the Mn are dilute) 
and vanishes as J ^ oo due to the contraction of the 
Bohr radius of the bound state. In this limit Tc depends 
crucially on the impurity band filling. For a full impu- 
rity band (one carrier per Mn) no low energy hopping 



processes are possible in a fully polarized ferromagnetic 
state; the ground state for a nearly filled impurity band 
is antiferromagentic or phase separated. The static mean 
field theory does not capture this physics at all, pre- 
dicting instead aTc ^ for all J. 

We now present a dynamical mean field theory which 
gives a reasonable account of the small and interme- 
diate J regime as well as the crossover to the 'impu- 
rity band' regime. It however does not adequately treat 
the band narrowing arising from extreme wavefunction 
localization so breaks down at some J >> Jc- We 
model the GaAs:Mn system as a lattice of sites, which 
are randomly nonmagnetic (with probability 1 — x) or 
magnetic (with probability x). Standard arguments p^ ] 
show that the relevant physics may then be determined 
from the local (momentum-integrated) Green function 

Groc(^) - b^Il^i^ - - ^pa)-\ Gicc is in 

general a matrix in spin and band (not shown) indices 
and depends on whether one is considering a magnetic 
(a) or non-magnetic (6) site. Being a local function, it 
is the solution of a local problem specified by the par- 
tition function Zioc — J dSe"^'"" with action Sioc = 

50a/3('r-T')cL(T)Ca/3(r') + JS-X;aa/3 cta{x)aaflCaf3{x) On 

the a (magnetic) site and Sioc = 5oa/3('^^'^')cL('r)ca/3(T') 
on the non-magnetic (6) site. The a-site mean field 
function can be written as ffg^^ — ao + aifh ■ Gap 
with m the magnetization direction and oi vanishing 
in the paramagnetic state. It is specified by the con- 
dition that the local Green function computed from Zioc: 
namely Sin Zioc/SgQ = (gg ~^)~^ is identical to the lo- 
cal green function computed by performing the momen- 
tum integral using the same self energy. The momen- 
tum integal requires an upper cutoff because the /2m 
dispersion given in Eq. |l| applies only near the band 
edges. We take the density of states as a semicircle, 
N{e) = b^J (0T'5(£-£p) = V4t2 - e727rt2 with param- 
eter t chosen to match the band-edge density of states, 
t = {AT:f/^ /2mb^ {t = 0.75eV and 1.5eV for the GaAs 
heavy and light hole bands). This choice of cutoff corre- 
sponds to a Bethe lattice in infinite dimensions; the cru- 
cial point is that it has the physically correct band edge 
density of states, this is the only aspect important for 
our work. Then gQ obeys the equation go('-^) — Sq{^) = 
u; + ^l- xt'^iig^iio) + JS ■ ac^p)-^) - (1 - x)egl{Loy^ 
where the angular brackets denote averages performed in 
the ensemble defined by the appropriate Zioc- 

The solution of the equation depends crucially on J/i, 
X and T. The inset to Fig. |^ shows the majority-spin 
density of states corresponding to the T = ferromag- 
netic state. For small J we see the expected shift pro- 
portional to xJ . For J > Jc = t an impurity band cen- 
tered at ~ — J and containing x states is seen to split off 
from the main band. The DMFT Jc is in good numerical 
agreement with the results of ||l^] ; this and the obviously 
correct qualitative behavior confirms its reliability in the 
experimentally relevant regime. 
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FIG. 1. Main panel: density of states at T > Tc for 
J = 0,0.5t,t and 2t. Inset: majority-spin density of states 
at T = for the same parameter values. 
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FIG. 2. Calculated transition temperature vs carrier con- 
centration for X = 0.05 and various J values as shown 

As the temperature is increased, the spins disorder 
and eventually the magnetic transition temperature is 
reached. Above this temperature, go is spin- independent. 
The main panel of Fig. shows the density of states for 
T > Tc. For xJ"^ /{t^ - J^) << 1 there is a small spin 
independent band offset of size xJ^ / {t^ — J^). For J > t 
an impurity band forms, corresponding to carriers locally 
parallel to Mn spins. 

The ferromagnetic transition temperature Tc may be 
obtained by linearizing the equation in the magnetic part 
of go, leading to an implicit equation for Tc- 

^ -2e{xjf/i 

^ ial x^PY (1 - t^lgl) - xJH^{b/i - PIgl) 

(4) 

where temperature is contained in the Matsubara sum 
over the frequency a;„ on which g^ depends. 



FIG. 3. Calculated Tc as a function of the local exchange 
coupling J for a fixed value (x — 0.05) of the Mn concentra- 
tion and for different values of the hole density per Mn ion 
{uh = 0.2, 0.5, 0.75, 0.85, 1.0 as shown). 

Figure. ^ shows the electron density dependence of 
the magnetic transition temperature for J = O.bt (less 
than the critical value for impurity band formation), 
J = t (the critical value for impurity band forma- 
tion) and J = \.bt and 2t (where the impurity band 
is well formed). The striking feature, evident in all 
three curves, is the non-monotonic behavior of the tran- 
sition temperature. This has different origins in different 
regimes. For x,P / (t^ — P) < 1 an analytic solution for 
Tc may be obtained. The details will be presented else- 
where; one result is that the density ti™^,'^ at which Tc is 
maximized is n""''' - ^(2 - .y3 + 2.P/t^ - J4/t4)3/2 _^ 
0{x,P/ {t^ - .P))) K. 0.04 - 0{J/t) . Thus in this limit 
the Tc maximum is a consequence of structure in the 
underlying electronic susceptibility. The precise posi- 
tion depends on the cutoff, but is very low. For J > t 
the physics is dominated by the spin-polarized impurity 
band. In this limit Tc is controlled by the delocalization 
energy in the impurity band, and is therefore maximized 
when the band is half filled. In a filled impurity band 
(n = x) no low energy hopping processes are allowed 
in a ferromagnetic state whereas in an antiferromagnetic 
state hopping is allowed with amplitude x^/^i/A where 
A ~ J is the gap between the impurity and conduction 
band. This physics implies that very near the filled impu- 
rity band limit, the ground state is antiferromagnetic. As 
A increases the window of antiferromagnetism decreases. 

Fig. ^ shows the magnetic transition temperature as a 
function of magnetic coupling J for different hole densi- 
ties ranging from rt/, = 0.1/Mn to 1/Mn ( in our con- 
ventions, n = xrifi) The collapse in Tc for the filled band 
is evident. The physically evident decrease of Tc at very 
large J due to decrease of impurity state Bohr radius is 
not captured by our model, so we expect that for J > 2t 
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FIG. 4. Variation of Tc with n at various Mn concentra- 
tion X and J — 1. 

our calculation overestimates Tc- 

Fig 4 shows the dependence of Tc on Mn concentration 
X for J = t, a, value of the order of the LSDA estimate. 
We see that simultaneous increases in the Mn concentra- 
tion (by say a factor of 2) and the density (by say a factor 
of 4) should increase by more than a factor of two. 

We compare our results to the predictions of other 
means of calculation. The 'mean field theory' ap- 
plied to our model predicts Tc = xn^^'^J^/t at all n,J. 
In the limit xJ^/ {t^ - J^) < 1 and < n< n"^"^"^ our 
analytic solution of the equations yields the mean field 
result but we find deviations either as J approaches 
or when n exceeds rv^'^^ or in the extremely low density 

limit n < {xJ"^ /tf'^^ (not visible in the plots shown here) 
where Tc ~ xn. An alternative approach to Tc involves 
spin- wave excitations |^J^. In classical high-spin mag- 
nets at T Tc spin waves are excited throughout the 
Brillouin zone and Tc occurs when the number of exci- 
tations (set by T divided by a typical magnon energy) 
is large enough. The present theory may be thought of 
as a calculation of a typical (i.e. averaged over the zone) 
magnon energy (which is itself determined by the changes 
in electronic energy due to spin disorder) on the assump- 
tion that the spin wave excitations have no particular 
spatial structure. In d = 3 for T near Tc this is cor- 
rect except for small amplitude critical fluctuations of no 
particular energetic significance. More importantly, our 
calculation provides detailed access to the experimentally 
relevant intermediate J,n,x regimes. 

We now briefly discuss numerical estimates of Tc for 
GaAs:Mn. The band theory estimates J = leV and 
t = Q.75eV (corresponding to the heavy hole mass) along 
with the typical density rih =0.1/Mn, implies a J/i ~ 1.3 
and a Tc ~ 80K for x = 0.05 whereas the light hole mass 
(t = l.beV) implies J/t ^ 0.75 and a Tc ^ 50K. These 
Tc values are for the single band model; the contributions 
of the two bands add, leading to a Tp'^v' « 130K. Al- 
though our theoretical estimates agree well with the ex- 



perimental Tc [|i],D, this agreement should not be taken 
too seriously in view of the simplifying approximations 
of our model. Interestingly, for the LSDA J the heavy 
hole band makes a higher contribution to Tc. Increase 
of n by about 50% will increase Tc by a similar amount 
for the heavy hole band but less for the light hole band. 
Increases in J if it can be managed will also increase Tc 
(although much less rapidly than the quadratic depen- 
dence predicted by the mean field theory). The most 
promising route to a higher temperature ferromagnet is 
predicted to be a simultaneous increase in cc to a value of 
order 0.1 and n to about 0.02 or about 0.2/Mn. 

In summary, we have presented a theory of the mag- 
netic semiconductors which can handle both the weak 
coupling limit (xJ less than Fermi energy Ep) and the 
intermediate coupling regime (J > t but not too large). 
This method correctly treats the physically relevant sit- 
uation in which the carriers are constrained to be locally 
parallel to the Mn spins and allows, for example, calcu- 
lation of the resistivity and optical conductivity. Discus- 
sion of these quantities will be given elsewhere. 
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